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Laboratory experiments measuring neutrino oscillations, indicate small mass differences between 
different mass eigenstates of neutrinos. The absolute mass scale is however not determined, with 
at present the strongest upper limits coming from astronomical observations rather than terrestrial 
experiments. The presence of massive neutrinos suppresses the growth of perturbations below a 
characteristic mass scale, thereby leading to a decreased abundance of collapsed dark matter halos. 
Here we show that this effect can significantly alter the predicted luminosity function (LF) of high 
redshift galaxies. In particular we demonstrate that a stringent constraint on the neutrino mass can 
be obtained using the well measured galaxy LF and our semi-analytic structure formation models. 
Combining the constraints from the Wilkinson Microwave Anisotropy Probe 7 year (WMAP7) data 
with the LF data at z ~ 4, we get a limit on the sum of the masses of 3 degenerate neutrinos 
£m„ < 0.52 eV at the 95 % CL. The additional constraints using the prior on Hubble constant 
strengthens this limit to Em„ < 0.29 eV at the 95 % CL. This neutrino mass limit is a factor ~ 4 
improvement compared to the constraint based on the WMAP7 data alone, and as stringent as 
known limits based on other astronomical observations. As different astronomical measurements 
may suffer from different set of biases, the method presented here provides a complementary probe 
of £m„. We suggest that repeating this exercise on well measured luminosity functions over different 
redshift ranges can provide independent and tighter constraints on Em„. 

PACS numbers: 98.80.Es, 14.60. Pq, 98.62.Ve, 95.80.+p 



I. INTRODUCTION 

A cosmic background of neutrinos is one of the key predictions of standard cosmology. Their predicted abundance is 
comparable to that of the relic photons. Thus if neutrinos have a mass they can contribute significantly to the matter 
density in the universe. Experiments which detect neutrino oscillations have measured small but non-zero differences 
between the mass eigenstates of neutrinos, with at least one of them having a mass larger than ~ 0.05 eV [l|, 0]. 
The absolute mass scale of neutrinos could be inferred from various /3-decay experiments [1, 0| ■ However at present 
stronger constraints on neutrino mass are obtained from cosmological observations. In particular observations related 
to anisotropies in the Cosmic Microwave Background Radiation (CMBR) and the growth of structure in the universe 
can play an important role (see 0, Q for reviews). For example the recent WMAP7 data itself has been used to set 
an upper limit on the of sum of neutrino masses, Em^ < 1.15 — 1.3 eV, for sudden reionization @, @] and Em„ < 1.7 
eV for a generalized reionization scenario [H . Thus one is dealing with the universe having matter density dominated 
by cold dark matter (CDM) with massive neutrino (Hot dark matter, HDM) providing sub-dominant contribution 
(i.e the neutrino density parameter fl u = Em„/93.14/i 2 eV < 0.014/i -2 , where h is the Hubble constant Ho in units of 
100 km s _1 Mpc -1 ; see chapter 2 of Ref. [H[). This is usually referred to as the mixed dark matter (MDM) scenario. 

Interestingly, the presence of even such a subdominant HDM component leads to the suppression of the growth of 
density perturbations below a scale known as the free streaming scale 0]. This suppression occurs because neutrinos 
escape (or free-stream out of) their own density perturbations below the free-streaming scale. Thus only the CDM 
component can lead to perturbation growth below this scale. This free streaming scale is time dependent and also 
depends on the neutrino mass. The suppression of the growth of density perturbations below the free streaming scale 
leads to a reduction in the matter power spectrum and a delay in the formation of structures in the universe. This 
in turn results in a reduced abundance of dark matter halos at any given epoch, above a characteristic mass scale. 
Thus observations related to large scale structure formation in the universe can be used to probe the absolute mass 
scale of neutrinos. It has been shown that sub eV constraints are obtained on neutrino mass by combining CMB data 
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with currently available data on galaxy surveys at low redshift counts of low redshift massi ve g alaxy clusters 
plj . inter-galactic medium (IGM) Lyman-a absorption power spectrum [l2|, EH and weak lensing 14 1. Unlike CMB 
observations all the astrophysical observations arc affected by different systematics related to baryonic physics. Thus 
it is important to explore whole range of observables (with different biases) to get realistic constraints on Em„. There 
is a growing wealth of observations on high redshift galaxies which may also provide independent and equally useful 
constraints on neutrino mass. Thus in this work we explore the possibility of using the luminosity functions (LF) of 
high redshift Lyman break galaxies (LBGs) for constraining Em„. 

The basic idea is as follows: The reduction in the matter power spectrum in models with Em„ > 0, compared to 
models where Em„ = implies a reduced abundance of galactic scale dark matter halos at high redshifts. In order 
to account for the observed LF of these galaxies (in number per unit volume per unit luminosity), the light to mass 
ratio of each galactic halo has to be systematically higher in the models with Xm„ > 0. However changing the light 
to mass ratio is degenerate with the unknown extinction correction one applies to the observed LF. Nevertheless, this 
degeneracy can be lifted if one has a feature in the LF at some characteristic mass scale, introduced by various feedback 
processes like the radiative feedback after reionization. In such cases the shape of the predicted galaxy luminosity 
function depends on the neutrino mass. We use this idea to constrain neutrino mass. In the next section we discuss 
structure formation models incorporating massive neutrinos. Section Mil describes our semi-analytic models for the 
UV luminosity functions. The effect of a non-zero neutrino mass on the high redshift LFs is studied in Section IIVI 
We present our limits on Em^ in Section [V] using a Markov Chain Monte Carlo analysis and conclude in Section PVTl 



II. STRUCTURE FORMATION IN MDM COSMOLOGY 



A crucial ingredient of any model of structure formation is cr(M, z), the rms density fluctuations on any mass scale 
M as a function of redshift z. This is given by 

a 2 (M, z) = a 2 (R, z) = J ^ k 2 P k (z)W 2 (k, R) (1) 

where R is the comoving radius of a sphere containing mass M, k is the comoving wave number, W(k,R) is the top 
hat window function in Fourier space and Pk{z) is the linear power spectrum of the density fluctuations at z. For a 
universe with massive neutrinos Pk(z) = T 2 (k, z, z.i;m v )Pk(zi). Here Pk(zi) is the initial power spectrum at Zi and 
the function T(k, z, Zi\ m u ) is the matter transfer function in a ACDM universe with massive neutrinos. 

The transfer function in mixed dark matter (MDM) cosmology with neutrinos has been studied extensively 0, EE 
EH- Eisenstein and Hu [l6| give a fitting formula of the form 

T(k, z, z,; m v ) = T master (k; m v )D{k, z, z { ; m„). (2) 

The function D{k, z, z%\ TnJ) is the scale dependent growth factor of CDM, baryon and neutrino perturbations in a 
ACDM universe in presence of free streaming neutrinos and T mas ter is a master-transfer function. The explicit forms 
of both these functions are given by Eisenstein and Hu [lH (see also [HJ). Their fit is optimized for a total number 
of three neutrinos which include massive as well as massless species. Note that the standard model of particle physics 
predicts the effective number of neutrinos {N v ) to be 3.04 with the 0.04 coming from incomplete neutrino freeze-out 
and finite temperature effects around e + e~ annihilation Q. We have checked that this transfer function is in good 
agreement with T computed numerically using the 'CAMB' code [l!| (Lewis, Challinor & Lasenby 2000). We will 
also fix N v — 3 and adopt this transfer function while calculating the luminosity functions below. 

Given the transfer function and cr(M, z) one can estimate the abundance of dark matter halos, for example using 
the Press-Schechter (PS) approach. This also requires one to specify 5 C , the linearly extrapolated critical density 
contrast needed for collapse of a spherical top hat over dense region. For a flat universe with only CDM, the critical 
density contrast of collapse is 1.686 [2lll22l]. The calculation of S c for the ACDM model without any massive neutrino 
species, has been done before by [23j, who found that 8 C for this model is very nearly the same as that for the pure 
CDM model. However, to our knowledge, a corresponding estimate of 5 C does not exist for models with massive 
neutrinos. As the abundance of halos is exponentially sensitive to the value of 5 C when using the PS formalism, it is 
important to estimate it even in MDM models. We do this below using the spherical model for nonlinear evolution 
adopting a flat ACDM universe with massive neutrinos. 



A. Spherical model for nonlinear collapse in MDM cosmology 

The spherical model or top hat model [24| , gives the nonlinear growth of a uniformly overdense spherical region in 
a smooth background of expanding universe. In spherical model we study the evolution of physical density contrast 
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S(r, z) in real space directly rather than the evolution of it's Fourier components S(k, z). For this purpose, we assume 
a spherical uniformly over dense region in the background expanding universe. Below the free-streaming scale of the 
neutrinos, only the CDM and baryons can cluster due to gravity, although both the cosmological constant A and 
neutrinos will contribute to the expansion. Note that the mass associated with the free streaming scale is much larger 
than the galactic scales we consider below. Thus the clustering mass is solely contributed by the baryons and the 
CDM. This is also borne out by recent simulations of structure formation in an MDM universe by Brandbyge et al 

The gravitational acceleration of the shell of radius r (with initial radius fj) in the spherical model is then given by 



r = — - 



GM(zi 



„2 3 



--nG(p v (z)-2p A )r. (3) 



Here M is the total mass which can cluster within the spherical region, and hence does not include any contribution 
from the neutrinos. Explicitly we have 

M = inrfipcDMizi) + Pb(zi)){l + S cb ( Zi ))/3 = 4tit?(1 - U)p m (z t )(l + S cb ( Zi ))/3. (4) 

Here pcdm{z), P b{z), p v {z) and P a(z) are respectively the background densities of CDM, baryons, neutrinos and 
cosmological constant at redshift z. We have also defined the neutrino fraction /„ = p^/pm, with p m = Pcdm + Pb + Pv 
is the total matter density. Further S cb = (SpcDM + &Pb) / {pcdm + Pb) is fractional overdensity in the CDM + baryon 
component within the spherical region. Note that the total density contrast S m = Sp m /p m = S c b(l — /„), since 
5pv = inside the spherical region. Even though neutrinos do not cluster below their free-streaming scale, their 
uniform density does lead to a deceleration of the shell, through the p v {z) term in Eq. ([3]). Moreover the cosmological 
constant leads to a positive acceleration of the shell proportional to {p A + 3p\), where the pressure p\ = —p A . At 
any redshift z and p A {z) = Pc Qa(0) and p v {z) = /0 C £7„(1 + z) 3 , where fii is the present density of component H' in 
units of the critical density p c . 

The Hubble expansion rate at any redshift is given by 

H( z ) = H(o)(a m (i + z) 3 + n A ) 1/2 (5) 

where H(0) is the present value of Hubble constant. Noting that 

dt da dz 

dt = -r~rdz = — TT . . .„ rr, (6) 

dadz H(z)(l + z) 3 w 

we can numerically solve Eqn. @ for a given set of initial conditions to obtain trajectory of the shell r(z), and the 
redshift of collapse z c . We assume that initially the density contrast 8 cb {zi) = Si is small enough that the overdense 
region is expanding along with the background. Thus the initial velocity of the shell at radius r is Uj = H{z i )r{z i ). 

Once we know the redshift of collapse z c , we can evolve linear equations of perturbations using the same initial 
conditions to calculate the linearly extrapolated critical density contrast, 5 C at z c . This linear evolution of 5 cb for 
scales much below the free streaming scale is governed by 



Jcb 



2HS cb - 4ttG(1 - f v )n mPc 5 cb = 0, (7) 



whose solution can only be determined numerically for a A+MDM universe. 

In the upper panel of Fig. [TJ we give 8 C for z c = as a function of fl A , as determined by the above procedure, 
assuming zero mass for the neutrinos. We confirm the weak dependence of S c on Q\ (with S c differing from 1.686 
only by 0.5% for Q A = 0.7), as also found by (23|. Note that even smaller changes in 5 C (~ 0.006%) would obtain in 
a ACDM model for the higher collapse redshifts (say z c = 4) relevant for our work. In the middle panel of Fig. [TJ we 
show the S c dependence on £m„, for z c = 4. We have adopted the cosmological parameters of a fiducial model with 
SIa = 0.735 and h = 0.71. The change in S c is very small even in the case of universe with massive neutrinos for the 
Y>m v range constrained by CMBR observations. For example, with z c = 4, S c decreases only to 1.68 for £m„ = 1 eV, 
or only a decrease of ~ 0.4% from the canonical value. The bottom panel gives S c as a function of z c for £m„ = 1 
eV and for the cosmological parameters as above. At z c ~ now S c decreases to 1.671 (or a fractional change of 
~ 0.9%), due to the additional effect of the ft a at such low redshifts. Based on the above discussion, it is a reasonable 
approximation to take the fiducial value of S c = 1.686, when computing the abundance of halos at high redshifts, even 
in the presence of massive neutrinos. 
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FIG. 1: Dependence of 5 C on various parameters. The top panel shows <5 C as a function of Qa for a ACDM universe without any 
massive neutrinos, with z c = 0. The middle panel shows S c as a function of £m„ for ACDM universe with massive neutrinos 
where z c = 4. The cosmological parameters are chosen according to our fiducial model (Oa = 0.735, h = 0.71) and by keeping 
Ocda/ + + = 0.265. In the bottom panel the dependence of S c on z c is plotted for Em„ = 1 eV and for the cosmological 
parameters as above. In all the panels the solid horizontal (solid dark) line at 1.686 corresponds to S c for a universe with only 
CDM in it. The red dotted line shows the variation in 5 C with the parameter of interest. 



B. The abundance of halos 



The abundance of dark matter halos as a function of mass and redshift plays a crucial role in understanding various 
aspects of galaxy formation. This halo mass function can be obtained through analytical formulations as well as N- 
body simulations. The first analytical form of halo mass function was given by Press and Schecter [20l |. Subsequently 
several alternative mass functions have been suggested which better fit the results of N-body simulations, a popular 
choice being the Sheth-Tormen (ST) mass function [27||- These mass functions assume the initial density field to be 
Gaussian random. And also that a region of size R (containing a mass M) collapses when its linearly extrapolated 
density contrast reaches a critical density S c , as predicted by the spherical model of nonlinear evolution. We found 
above that S c for the CDM+baryon component is only slightly altered from the canonical value of 1.686, by the 
presence of massive neutrinos. 

Brandbyge et al [25| have shown that the ST mass function provides an excellent fit for mass functions of dark 
matter halos obtained from N-body simulations including massive neutrinos. This has also been confirmed recently 
by Marulli et al [26| . In the ST formula the comoving number density of halos with mass between M and M+dM at 
any redshift z is given by 

where from the numerical fitting A — 0.3222, a = 0.707, p = 0.3 and S c = 1.686. Brandbyge et al [25| also find 
that one needs to use p = p c b = (pcdm + Pb) and calculate the mass as M = M c b = ^r 3 {pcDM + Pb)/3 in using 
ST formula to fit the halo abundance. However they calculate a using the total matter power spectrum (25j (also 
Brandbyge; private communication), although naively one may have expected that a should be calculated using the 
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FIG. 2: The halo abundance Nst for ACDM and AMDM models. The top panels show Nst as a function of halo mass for 
£m„ = (solid black lines), Em„ = 0.5 eV (long dashed red lines) and Sm„ = 1 eV (dash-dotted blue lines). The bottom 
panels show the ratio of Nst for different values of Sm„ to that obtained taking Em„ = 0. The left panels are for z = 0, while 
the right panels are for z — 4. 



power spectrum of CDM+baryons 1 . In our work we follow the prescription of [25[ in calculating the halo abundance, 
as this seems to fit the simulated data well. 

We illustrate the effect of massive neutrinos on the abundance of dark matter halos in Fig. [2] where we plot in the 
upper panels, Ngx, using the prescription of (25|, for two redshifts z = and 2 = 4. The solid black line in Fig. Ogives 
Nst, for a fiducial ACDM model with Hm v = 0. As the fiducial ACDM model, we adopt a set cosmological parameters 
consistent with WMAP-CMBR 7 year data [H|. Thus we adopt h = 0.71, Q b = 0.0448 and n CDM = 0.220. The 
initial power spectrum is specified by the best fit values of scalar spectral index n s = 0.953 and curvature fluctuation 
amplitude = 2.43 x 10~ 9 at a pivot scale fco = 0.002Afpc _1 (or a as = 0.801). We then examine models with 
non-zero ~Em u by varying Qcdm and satisfying the constraint that Qcdm + £l v = 0.220. The long dashed (red) 
line gives Nst for the case when Em^ = 0.5 eV, while the dash-dotted (blue) line is for the case T,m l , = 1 eV. The 
bottom panels give the ratio of Nst between the non-zero neutrino mass cases with respect to the fiducial model, as 
a function of M. 

From this figure we can see that the number of collapsed halos decreases with increase in neutrino mass for halo 
mass scales above a characteristic mass scale M c , while there is a slight increase in abundance of halos below this 
mass scale. The characteristic mass M c ~ lO 12 Af at z = while M c ~ 1.6 x 1O 6 M at z = 4 for m v = 1 eV. 



Note that, for length scales much smaller than the free-streaming scale, since S m = S cb (l — /„), a calculated from the total matter power 
spectrum is ss o" c (,(l — f v ), where is that calculated from the power spectrum of the CDM+baryons. 
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The decrease in abundance of large mass halos arises as a result of the suppression of growth due to the presence of 
free-streaming neutrinos. In our model with Em„ = 1.0 eV, we find that at z = 0, the abundance of galaxy clusters 
of mass M ~ 10 14 M Q decreases by a factor ~ 3.2, while for bigger clusters with M ~ 1O 15 M0, the decrease in 
abundance is by a factor ~ 25. Similarly at z = 4, the abundance of galactic scale halos of mass 10 11 M Q and 10 12 M© 
are suppressed respectively by factors ~ 6.5 and ~ 30, in models with Hm v = leV compared to the zero neutrino 
mass case. Thus strong constraints on the neutrino mass can be obtained from measuring the abundance of galaxy 
clusters at low redshifts or equally, the abundance of galaxies at high redshifts as probed by their observed luminosity 
functions. 



III. HIGH REDSHIFT GALAXY LUMINOSITY FUNCTIONS 



In the previous section we have described how the presence of massive neutrinos suppresses the high redshift galactic 
halo formation. In this section, using semi-analytical models, we investigate the effect of massive neutrinos on high 
redshift galaxy UV LF. We model high redshift galaxy LF using the semi-analytical treatment of Samui, Srianand & 
Subramanian (2007) (3l| (hereafter SSS07) which were successful in explaining the observed luminosity functions at 
3 < z < 10. (see also Samui, Subramanian & Srianand (2009) [3(|; hereafter SSS09). Here, we use the same approach 
in order to constrain neutrino masses. We briefly describe the modeling here and the interested reader may refer to 
SSS07 and SSS09 for more details. 

The star formation rate of an individual dark matter halo of mass M collapsed at redshift z c and observed at 
redshift z is modeled by (see HH, [33l ]) 



M S f(M,z,z c 



1 L rn _ 



*(*) ~ t(Zc) 

K dyn(z c ) 



X exp 



t(z) - t{z c ) 



^dyn(z c ) 



(9) 



where /» is the fraction of the total baryonic mass that is converted into stars over the entire lifetime of the galaxy. 
Here td yn (z c ) is the dynamical time scale of a halo collapsing at z c (Eq. (3) of SSS07) and k is a parameter which 
governs the duration of the star formation activity in the halo which we take here to be unity. Further, t(z) is the 
age of the universe at redshift z; thus t(z) — t(z c ) is the age of the galaxy at z. We assume that stars are formed 
with a normal Salpeter IMF in the mass range from 1 — 100 M Q . The population synthesis code Starburst99 [34[ 
is used to obtain the luminosity of a galaxy undergoing a burst of star formation, as a function of time for a given 
rest wavelength of 1500 A. The assumed star formation rate of a galaxy in Eq. (|9]) is then convolved with this burst 
luminosity to get the evolution of luminosity (£1500) of an individual star forming galactic halo (See Eq. (6) and 
figure 1 of SSS07). Only a fraction (I/77) of the total light produced by the stars comes out of the galaxy due to the 
absorption by dust. We convert this luminosity (L = Li^qo/tj) to a standard absolute AB magnitude M A b using the 
equation given by (35| . The luminosity function <&(Mab, z) at any redshift z is then given by (3l| 



$(M AB ,z) dM AB = / d 



dN ST (M, z c ) 
dz c 



dM dL 



1500 



dL 



1500 



dM 



dM AI 



AB 



(10) 



where dNsr{M, z c )/dz c = Nst{M, z c )dt/dz c , and Nst(M, z c )dM is the formation rate of objects in the mass range 
(M, M + dM) at redshift z c . We model this formation rate as the time derivative of ST mass function as (i) they are 
found to be the best in reproducing the observed LF of high- z LBGs (sec SSS09), and (ii) as explained earlier gives 
good fit to the abundance of dark matter halos in N-body simulations incorporating massive neutrinos [25| . 

Star formation in a given halo also depends on the cooling efficiency of the gas and various feedback processes. We 
assume that gas in halos with virial temperatures (T V i r ) in excess of 10 4 K can cool (due to recombination line cooling 
from hydrogen and helium) and collapse to form stars. However the ionization of the IGM by UV photons increases 
the temperature of the gas thereby increasing the Jean's mass for collapse. Thus in ionized regions, we incorporate 
this feedback by a complete suppression of star formation for halos with circular velocity v c ^ 35 km s _1 and no 
suppression with v c ^ V u = 95 km s _1 13611 . For intermediate circular velocities, a linear fit from 1 to is adopted 
as the suppression factor ([IH; see also }37t l38l]. SSS07). In SSS07 we found that this feedback mechanism naturally 
leads to the observed flattening of the LF at the low luminosity end 2 . Note that implementing radiative feedback in 



2 In principle, one could also consider star formation in smaller mass halos with T v i r > 300 K, where the cooling is due to the H2 
molecular line emission [39I ]. While star formation in such halos could be important for rcionization and IGM metal enrichment, in the 
post reionization era, it is suppressed by radiative feedback. 



7 



z 


Em„ = eV 


£m„ = 0.5 eV 


Em„ = 1 eV 


f*h x 2 xl 


f*/v x 2 x 2 


Uh x 2 x 2 


3.0 
4.0 
5.0 
6.0 


0.040 15.14 1.89 
0.039 16.09 1.46 
0.030 45.71 5.07 
0.041 4.34 0.72 


0.052 19.42 2.16 
0.067 19.49 1.77 
0.054 44.87 4.99 
0.082 5.64 0.94 


0.086 31.39 3.92 
0.152 40.08 3.46 
0.179 12.91 1.43 
0.250 4.97 0.83 



TABLE I: The best fit values f*/r), corresponding \ 2 an d reduced xl f° r our models at various redshifts for three different 
values of Em„. 



a model requires a knowledge of the reionization epoch. A number of observations suggest that reionization is nearly 
complete by a redshift z re > 6 (e.g. 0, EH)- We will show below that the strongest constraints on neutrino masses 
are obtained from the luminosity function at z = 4. As this redshift is much lower than z re , our limits on the neutrino 
mass then turns out to be insensitive to the exact z re as long as it is greater than 6. 

In our models, we also incorporate the possible Active Galactic Nuclei (AGN) feedback that suppresses star for- 
mation in the high mass, by multiplying the star formation rate by a factor [1 + (M /Magn) 3 ]' 1 , as in SSS07. This 
sharply decreases the star formation activity in high mass halos above a characteristic mass scale Magn, which is 
believed to be ~ 10 12 M Q (see [13, S3)- The determination of Magn is done by fitting luminosity functions at a 
fiducial redshift and will be discussed below. In what follows, we compute the high redshift luminosity functions 
of galaxies in the presence of massive neutrinos using the above semi-analytical prescription and comparing it with 
observations to put limits on neutrino mass. 

IV. EFFECT OF NEUTRINO MASS ON HIGH REDSHIFT LF 

We now illustrate the effect of massive neutrinos on high redshift galaxy UV LF. We show in Fig. [3J our model 
predictions of luminosity functions at different redshifts, z = 3 — 6, along with the observational data. The redshifts 
arc indicated at the top of each panel. The observed data points are taken from Reddy et. al (2007) [3 for z = 3, 
and Bouwens et. al. (2007) [3 for z = 4, 5 and 6, and corrected to the fiducial WMAP7 cosmology described earlier. 
As in SSS09, we take into account cosmic variance (13; 13) an d add an uncertainty of 14% to the Poisson error in 
quadrature for redshift 4 — 6 data. 

The solid line shows the predicted best fit luminosity functions at various redshifts for the fiducial cosmology and 
with Y>m v = 0. In order to fit the observed data points we have adjusted the free parameter f*/r] in our models, using 
X 2 minimization. We use all the points excluding the last two data points for z = 4, 5 and the last point for z — 6, 
in the low luminosity end while fitting the LF, as these points are affected by the incompleteness of the survey (see 
3)- For this illustrative study, we also fix V u = 95 km s _1 and adopt Magn = 1.8 x 10 12 M Q that best fits the LF 
at . z = 4 with Em„ = 0. 

Our model with £m„ = (solid curve) reproduces the observed LF very well at all the redshifts except z = 5 
(where SSS09 also earlier pointed out one discrepant point). The best fit values of f*/r) at different redshifts and 
the corresponding \ 2 an d reduced xli are tabulated in Table [I] The required values of f*/r), and reduced \li at 
different z arc nearly the same and similar to what we obtained in SSS09 (which used slightly different cosmological 
parameters). The flattening of the predicted LF as seen in Fig. [3] at the faint end is due to the radiative feedback. 

The thick dashed-dotted (blue) curves show the predicted LF if we use the same /*/?? for a MDM model with 
YiTTiy = 1.0 eV. It is clear that there is an order of magnitude suppression in the number density of galaxies at a 
given luminosity, which increases with increasing redshift. This is because the presence of neutrinos suppresses the 
formation rate of halos at the mass and redshift scales of our interest. We can make our model predictions match 
with the observed data by increasing f*/r] (i.c shifting this curve along the luminosity axis). These best fit values of 
f*/i) at different redshifts for the case £m„ = 1 eV, and the corresponding reduced xli are also tabulated in Table fl] 
Firstly, we can see from Table |H that at each z the value of f*/r/ (or light to mass ratio) needed to fit the observed 
luminosity function increases with Y,m v . For example, at z = 6 one needs ~ 6 times more baryons to convert into 
stars for a model with £m„ = 1 eV compared to the zero neutrino mass case. Furthermore, for Em,, = 1 eV, the 
value of f*/r] required to fit the UV LFs also increases significantly with z, in contrast to the zero neutrino mass case. 
Therefore, any independent constraint on f*/rj, especially at high z, could lead to useful constraints on the neutrino 
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FIG. 3: UV LF of LBGs at redshifts 3, 4, 5 and 6. The solid (black) line shows the predicted best fit LF for our model with 
Em„ = 0. The thick dashed-dotted (blue) curve shows our model predictions with Sm„ = I eV and using the same best fit 
f*/rj as the Em„ = eV model. This is to illustrate the suppression due to massive neutrinos. The thin dashed-dotted (blue) 
curves are best fits for the models with £m„ = 1 eV. The data points (filled and open triangles) for z — 3 are taken from [4l| 
and for z = 4 — 6 are from [4^] . 



mass. 

More importantly, these best fit luminosity functions, obtained with the new f*/rj (thin dashed-dotted curves), have 
a very different shape compared to the zero neutrino mass case. In particular, the predicted luminosity functions in 
models with £to„ = 1 eV, are suppressed at the low luminosity end compared to the zero mass case. This is basically 
because increasing f*/r] increases the light to mass ratio, which brings even small mass galaxies whose star formation 
has been suppressed due to radiative feedback, into the observable luminosity range. Therefore strong constraints on 
the neutrino mass can in principle be obtained by comparing the shape of the predicted luminosity functions with 
observations, independent of the free parameter /*/?7, provided the observations cover a wide luminosity range and in 
particular extend to the very faint end of the LF. 

This requirement is at present best satisfied at z = 4. This is because, the observed LF for z = 4 is well defined 
over a wide range of luminosity thanks to the Hubble Ultra deep field (HUDF) and the Great Observatories Origins 
Deep Survey (GOODS) (@, SE|). This redshift is also well below the redshift of reionization 0, S3]; so radiative 
feedback can be applied without any ambiguity. Therefore for demonstrating that quantitative limits on £m„ can be 
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Parameter WMAP7 



WMAP7 WMAP7 WMAP7 WMAP7 WMAP7 WMAP7 WMAP7 WMAP7 
+LF +HST +LF (m„) +LF +LF +LF(FB) LF(FB) 

+HST (m„) +HST +HST 

K) 



w 2 n b h 2 
ion DM h 2 



f*/v 



990+0. 060 
z.zzo_ 05g 

1 1 7 9+0.070 
1 - ±,z -0.068 

0.087±°:°S 
0.962t° ;™ 

o.7i7±S:8S 

66.1±li 



< 0.091 



9 99Q+O.O55 
^•^^»_0.056 

1 990+0. 065 
±.zzo_ 064 

090+ - 013 

u.uyu_ 014 
Qfi1+ 0015 

U-»Ol_0 .014 

820+ 051 
u.ozu_ 052 



H2.8 
-2.8 



65.5 



ns5 +0 007 

< 0.048 



9 9 R q+0.050 
z.zuo_ 056 

1 107 +o.05i 

■ L - lu '-0.050 

091 + 015 

n 070+0.013 

U-9(O_ 013 

7Sfi+ 049 



< 0.052 



9 9PQ+0.054 

Z.ZO»_q q 54 

i-i69l°;°58 

0.095±g : g« 
0.969t°:°lS 
0-821±g;gg 

036+ 007 
< 0.028 



9 9z ir+0.060 



1.177 



+0.071 
-0.069 



0.0861°;°!* 
0.9621°:°^ 
0.717t°:°£ 



66.0 



,+4.2 
4.0 



9 99Q+O.O55 
z.zz»_ 056 

1 9OC+0.065 
±.ZOO_q Qg 4 

0.090±g;SiI 
0.961±g;g^ 
0.824t°;°£ 



65.3 



H2.9 
-2.8 



n 0S4+ 005 



9 9 r 7 +0.053 
Z.ZUf _o.053 



1.168 



+0.056 
-0.058 



0.095±g;gi| 
0.970±g;gig 
0.821±g : gg 



68.7 



+2.2 
2.2 



n35 +0 007 



9 9or+0.056 
z.zoo_ 055 

i 9 na +0.060 
1.ZU8_ 05g 

0.090ig;gl| 
0.963lg;gl* 
0.808±g : gH 



2 266+ 051 

z.zuu_ 052 
l.lOO_ Q45 

093+ 014 
u.uyo_ 015 

? +0.013 

u.y/i_o.oi2 
n si q+0.030 

U - OIJ -0 .029 



66.3i5:| 
n q 7 +o.oo5 

< 0.042 



69.2 



+2.2 
-2.2 



r, o 7 +0.005 

u.uoi _ 005 
< 0.026 



Sm„(eV) < 1.0 



< 0.55 



< 0.54 



< 0.31 



< U 



< 0.52 



< 0.29 



< 0, 



< 0.28 



TABLE II: Results of our MCMC analysis to constrain Sm„. The first column lists the set of parameters that are obtained 
from our MCMC analysis. The last two rows give limits on /„ or Em„, at the 95% CL. For cases where we constrain f v , 
the neutrino mass is given by Em„ = 93.14 f v Qd Aih 2 , where for flnAih 2 we use the mean value given in row 3. For other 
parameters their mean values and la range are given. 



obtained with well defined LF, we will concentrate below on the LF of LBGs at z = 4. 

At this redshift we already see from TablcUthat, while \ 2 = 16.09 for Ym v — eV model, it increases by A\ 2 ~ 3.4 
for Em^ = 0.5 eV and by A\ 2 ~ 24 for £m„ = 1 eV models. This suggests that Em^ ~ 1 cV is strongly disfavored 
at a 5(7 level, and a typical 2a upper limit is close to Sm„ ~ 0.5 eV. While this is encouraging, we have been using a 
fixed set of cosmological parameters and therefore one needs to check if such a conclusion also follows when we vary 
these parameters. We examine this issue further below. 



V. LIMITS ON NEUTRINO MASS 



In the previous section we adopted a fiducial cosmology to examine the effect of a non-zero neutrino mass on the 
LF of LBGs. In order to obtain quantitative upper limits by exploring the full range of cosmological parameters 
consistent with both the WMAP7 data and the observed LF of LBGs, we have performed a Markov Chain Monte 
Carlo (MCMC) analysis using the publically available CosmoMC code [46[ (we use the CosmoMC May 2010 version 
and the WMAP likelyhood code version 4.1). The default version of CosmoMC constrains the neutrino fraction in 
dark matter defined as /„ = £L v /£Idm where Qdm = &CDM + &v So we constrain f v and obtain neutrino mass as 
Yimy = 93.14/„f2rjM^ 2 , where we use the mean value of SIdm from our MCMC analysis. We will also carry out the 
MCMC analysis by giving a prior in terms of Em„ to directly constrain its value. In addition to /„ (or Em„), we 
explored the usual 7 dimensional parameter space (flbh 2 , fl c h 2 ,8 s , r, n s , A s ,Asz)- 

In order to combine the constraints from the UV LF data with other constraints (like the WMAP7 data), we 
have added a new module to CosmoMC. This module computes the likelihood for LF of LBGs by comparing the 
theoretically predicted LF for any set of cosmological parameters with observed LF. The observed data points at 
z = 4 by [45|, are for a cosmology with f2 m = 0.3 and Q\ = 0.7 and h = 0.7. Therefore for each set of cosmological 
parameters in the chain, we first correct the observed data to that cosmology, by modifying the volume and distance 
scales appropriately. Then for each member of the chain we compute the theoretical LF, compare it with the observed 
data, and minimize the \ 2 by varying f*/r]. This minimum x 2 is used for calculating the LF likelihood. 

We give the constraints on Hm v obtained from the MCMC analysis in Table [TT| and Fig [5J The convergence of 
MCMC chains were diagnosed using the usual Gelman-Rubin statistic and an R — 1 < 0.01 was achieved in general. 
First we have examined the constraints from WMAP7 data alone @, S3- In this cas e we obtain an upper limit 
f v < 0.091 at the 95% CL, which when converted to a mass limit as described above, gives T,m u < 1.0 eV. The other 
cosmological parameters, given in column 2 of Table [Til are almost identical to that obtained by Q for their case of 
sudden reionization. Note that constraining Em„ directly gives a very similar limit Em„ < 1.08 eV at 95% CL with 
almost the same values for the other cosmological parameters (see column 6 of Table [TTJ) . 

Combining the UV luminosity function data at z — 4 with the WMAP7 data gives a significantly lower limit on 
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WMAP7+HST 



WMAP7 + HST + LF. 





(a)Best fit Cls 



(b)Best fit LF at z=4 



FIG. 4: The fits to the CMBR TT power spectrum (Cls) (left panel) and the LF at z = 4 (right panel) using the best fitted 
values of cosmological parameters obtained from our MCMC analysis. Solid blue lines are obtained for cosmological parameters 
corresponding to WMAP7+HST+LF data. The red dashed line in the left panel gives for comparison, Cls for best fitted 
cosmological parameters based on only WMAP7+HST data. 



£m„. We get f v < 0.048 corresponding to £m„ < 0.55 cV at the 95% CL. Directly constraining the neutrino mass 
also leads to a very similar limit, Y>m u < 0.52 eV at the 95% CL. Thus the limit on the neutrino mass is decreased 
by a factor of ~ 2 by the addition of the constraints from the z = 4 UV LF. The cosmological parameters for these 
two cases, given respectively in column 3 and 7 of Table HJ are very similar. These parameters are also similar to the 
parameters obtained from WMAP7 alone (column 2 of Table ILl]) . except for erg, which is increased on including the 
constraints from the LF data. Note that these limits on Sm„, obtained by varying all the cosmological parameters, 
bear out the naive expectation from the analysis of the Section IIV1 where we fixed the cosmology. 

Next we examined the effect of adding the constraints from the H determination of the HST SHOES (Supernova 
HO for the Equation of State) program (which we refer to as HST) ([5(3]; see also [H]). First WMAP7+HST data 
alone leads to a tighter neutrino mass limit of /„ < 0.052 (or £m„ < 0.54) at 95% CL, consistent with the earlier 
results of [4jj. On adding in the constraint from UV luminosity function data at z = 4 this limit is also further 
decreased. We obtain f v < 0.028 corresponding to Em y < 0.31 eV. (The corresponding cosmological parameters are 
given in column 4 and 5 of Table |H] respectively) . Again directly constraining the neutrino mass gives a very similar 
limit, £m„ < 0.29 eV at the 95% CL. Therefore, we see that adding in the LF data further decreases the limit on the 
neutrino mass by another factor of ~ 2. This limit, Sm^ < 0.29 eV is almost a factor ~ 4 improvement compared 
to the limit obtained using the WMAP7 data alone. The other cosmological parameters for this case are given in 
column 8 of Table [TT] and Fig 2] shows the corresponding predicted TT power spectrum of CMBR (Cls) and the LF 
at z = 4 overplotted against the data. In Fig [4] we also plot the Cls corresponding to the parameters which do not 
include the LF constraint (column 4 of Table HI)) . which shows that Cls for these two cases are almost identical. 

All the above neutrino mass limits along with the other cosmological parameters are summarized in Table [TTJ The 
corresponding ID marginalized distribution for /„, from various MCMC analysis, is shown in the top left panel of 
Fig. O The top right 2 panels of Fig. [5] shows 68% and 95% marginalized distributions for Hq and VtuMh 2 against 
f v , the massive neutrino fraction of the dark matter. The results of directly constraining Y,m u is shown in the 
middle panels of Fig. [5] We see also from these figures adding the constraint from the z = 4 UV luminosity function 
significantly improves the constraint on neutrino masses. 

In the above analysis we fixed the feedback parameters to Magn = 1.8 x 10 12 M q and V u = 95 km s _1 , as described 
earlier. It is of interest to test the sensitivity of our neutrino mass limits to changes in these values. One way of 
doing this would be to vary Magn-, Vu and f*/rj as free parameters to obtain the best fit LF and its corresponding 
likelihood for each step of the MCMC analysis. However, implementing the above procedure in CosmoMC (where 
the cosmology is also varied) is computationally expensive, and is outside the scope of the present work. Note 
however that the feedback parameters are expected to depend on physical nature of the feedback mechanisms and 
not on cosmology. Therefore, as an alternative to varying the feedback parameters within CosmoMC, we adopt our 
fiducial cosmology and carry out the following procedure. We vary Magn, V u and f*/rj in the parameter ranges 
lO n M < Magn < 5 x 10 12 M q , 85 kms" 1 <V U < 120 kms -1 and 0.005 < f*/r) < 1, to obtain the best fit LF and 



11 




FIG. 5: Various ID and 2D marginalized distributions from our analysis. The top panel corresponds to our fiducial LF 
model with feedback parameters Magn = 1.8 x 1O 12 M0 and V u = 95 km s _1 . The middle panel shows plots corresponding 
to constraining m v directly with same feedback parameters as above. The bottom panel is for the LF model with feedback 
parameters Magn = 1.5 x 1O 12 A40 and V u = 105 km s _1 (model FB). In each row the left panel gives the marginalized ID 
distribution for /„ or m„. The vertical dashed lines correspond to 95% confidence levels. The other two panels in each row, 
show the regions of 68% (dark color) and 95% (light color) confidence levels for Ho and QdmIi 2 against /„ or m„. The various 
contours corresponds to constraints obtained using WMAP7-only (red), WMAP+LF (blue) and WMAP+HST+LF (green) 
data. 



its corresponding \ 2 as a function of Sm v . This best fit \ 2 as a function of Sm„, obtained after varying all 3 free 
parameters given above, is shown in Fig. [S] It is clear from Fig. [S]that even after varying the feedback parameters, 
the global minimum value of x 2 obtains for the zero neutrino mass case. This global minimum obtains for V u = 105 
km s~ x and Magn = 1-5 x 10 12 Mq. The typical 2a limit is Em„ = 0.35 eV, even in this case where all the feedback 
parameters are varied. This excercise suggests that the uncertainity in the feedback parameters may not have a strong 
influence on the the upper limit on Em„. 

We have repeated our MCMC analysis adopting these best fit values of Magn a nd V u (which we label as model 
FB). The results of the analysis, adopting the feedback parameters of model FB, are presented in the last 2 columns of 
Table |ll]and in the bottom panels of Fig. [5] Combining the WMAP7 data and the LF data now leads to a lower limit 
on the neutrino fraction /„ < 0.042 and a corresponding Ym v < 0.48 eV. On adding also the constraints from the Hq 
determination of the HST SHOES program (together with WMAP7 and UV LF), we get /„ < 0.026 corresponding 
to Em„ < 0.28 eV at the 95% CL. We note that these limits on Em,, are almost identical to that obtained using the 
fiducial feedback parameters. 

In passing we note that we can also get quantitative measurements of the astrophysical parameter f*/r) at z = 4 
from our MCMC analysis. These values are given in Table [n] while in Fig [7] we plot the ID marginalized probability 
distributions (PDF) of /*/?? obtained from the MCMC analysis for a number of models. The red curves are for the 
fiducial feedback parameters (Magn = 1-8 x 10 12 M Q and V u — 95 km s _1 ) and blue curves are for the model FB 
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FIG. 6: The Sm„ vs \ 2 curve for our fiducial cosmology with massive neutrinos. Here we varied Magn, V u and f*/rj as free 
parameters in our model for LF. The solid black curve gives best fit x 2 as a function of m v after marginalizing over all the 
above three parameters. The number of degrees of freedom is 9. 




FIG. 7: The ID marginalized distributions of /*/?? from our analysis using LF at z = 4. The solid and dotted (red) curves 
respectively gives the PDF of j*jr\ when one uses WMAP7+LF and WMAP7+HST+LF data with the feedback parameters 
as in the fiducial LF model. The blue dashed and dashed-dotted curves give the PDF of f*/rj when WMAP7+LF and 
WMAP7+HST+LF data are used, and for the LF model with feedback parameters FB. 



(Magn = 1-5 x 10 12 Mq and V u = 105 km s _1 ). In particular the solid (red) and dashed (blue) curves give the PDF 
of f*/i] when WMAP7 and LF data are combined. The dotted (red) and dashed-dotted (blue) curves are obtained 
when constraints from WMAP7, HST prior on Ho and LF are combined. From this figure and Table [Til we see that 
fsf/ri is constrained in a narrow range of values, with a mean f*/rj ~ 0.034 — 0.037. If one adopts a value of r\ ~ 4 at 
z = 4 [El]], then this corresponds to an /* ~ 0.14 — 0.15. Thus about 15% of baryons need to be converted to stars 
over several dynamical timescales, to explain the observed UV LF at z = 4. This conclusion obtains regardless of the 
feedback parameters and although we have now explored the full range of cosmological parameters consistent with 
CMBR, HST and the LF data. 

The results of our MCMC analysis in this section shows therefore that fitting the UV LF of high z galaxies can 
significantly strengthen constraints on the neutrino mass. 
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VI. CONCLUSIONS 

We have proposed here a novel probe of neutrino masses which uses the high redshift UV LF of Lyman break 
galaxies. In particular, our models constructed in the framework of MDM cosmology, show that the observed shape 
of the UV LF of high redshift LBGs can be used to constrain the mass of the neutrinos (Em,,). 

In the presence of massive neutrinos, the matter power spectrum is suppressed below the free-streaming scale. 
This in turn results in a suppression of the abundance of collapsed halos capable of hosting galaxies at high redshift, 
compared to a model where Eto„ = 0. We have studied how this affects the UV LF of LBGs from z = 3 — 6. The UV 
luminosity of a galaxy in our models depends on the parameter f*/r], where /» is the total fraction of the baryons 
converted into stars and 1 jr\ is the fraction of the total light which comes out of the galaxy after taking into account 
dust extinction. For a given f*/r], we find that the number density of galaxies at a given luminosity is suppressed for 
models with a non-zero neutrino mass compared to the zero neutrino mass case. For example, if Sm„ = 1 eV this 
suppression factor is about an order of magnitude at all redshifts, and increases with redshift. In order to best fit the 
observed UV LF for the models with non zero neutrino mass, the light to mass ratio of each galaxy, governed by the 
parameter f*/r] has to be correspondingly larger at any redshift. At z = 6 for example, we find that f*/rj has to be 
6 times larger for a model with Em„ = 1 eV compare to the zero neutrino mass case. Moreover, unlike the Em„ = 
case, the best fit f^/rj also increases with z, for models with a non zero Em„. Thus independent constraints on /*/ry, 
especially at high z, could lead to useful constraints on the neutrino mass. 

More importantly, these best fit luminosity functions with a non-zero Yim v , obtained by increasing the light to mass 
ratio of each galaxy (or f*/rf) brings even small mass galaxies, which suffer from various feedback suppression effects, 
into the observable range. We have shown here that this results in a sensitivity of the LF shape to the neutrino mass. 
The well measured UV LF of LBGs at z ~ 4 is best suited for obtaining quantitative constraints as the data extends 
over the widest range of luminosities and to much fainter levels than at other redshifts. To obtain quantitative upper 
limits on Em,,, that also explores the full ran ge o f cosmological parameters, we have carried out an MCMC analysis 
using the publically available CosmoMC code J46[, after adding also a module which calculates the LF likelihood. The 
neutrino mass limits obtained here and other cosmological parameters for these models are summarized in Table [TTJ 
The corresponding ID and 2D marginalized distribution for /„ (or Em„), from various MCMC analysis, are shown 
in Fig. 

First, combining the constraints from the WMAP7 data and the UV LF of LBGs at z ~ 4, wc obtain a constraint 
on the neutrino fraction /„ < 0.048 and a corresponding limit on sum of neutrino masses, Em v < 0.55 eV at the 95 
% CL. Directly constraining Eto„ in the MCMC analysis leads a similar upper limit Em„ < 0.52 eV at the 95 % 
CL. Thus the neutrino mass limit is not greatly sensitive to whether one specifies a prior on /„ or Em„. We have 
also tested the sensitivity of these limits to the feedback parameters. Adopting a different set of feedback parameters 
(model FB) we get f v < 0.042 or Em„ < 0.48 eV at the 95 % CL. Thus the uncertainties related to halo mass 
(or circular velocity) range over which suppression of star formation due to feedback takes place and exact value of 
Magn, do not introduce significant uncertainty in the upper limit on Em„. Our constraints on Em„, obtained by 
combining the z = 4 LF and the WMAP7 data, are a factor of 2 more stringent than constraints that would obtain 
by using the WMAP7 data by itself. 

We have also examined the effect of adding the constraints from the HST prior on Hq. Consistent with the 
earlier work of [491 ] . we find that WMAP7+HST data alone leads to a tighter neutrino mass limit of /„ < 0.052 (or 
Etoj, < 0.54) at 95% CL. This limit is also further decreased to /„ < 0.028 or Em„ < 0.31 eV on adding in the 
constraints from the z = 4 LF data. Directly constraining the neutrino mass gives in this case Em„ < 0.29 eV at the 
95% CL. Adopting model FB for the feedback parameters also gives a similar limit /„ < 0.026 or Eto„ < 0.28 eV at 
the 95 % CL. We note that these neutrino mass limits are almost a factor ~ 4 improvement compared to the limit 
obtained using the WMAP7 data alone. 

A summary of the current cosmological and astr ophysic al constraints on neutrino mass can be found in [52| . Some 
of the specific results on E77i„ limits are given by fiol - flil . I53T456T ] . The constraints on Ymi v obtained here adding in 
the z ~ 4 UV LF data to the WMAP7 and HST data, are comparable (or better in several cases) to the above limits. 
Our work is mainly a demonstrative first step, where we have suggested the utility of the LF of high redshift galaxies 
to constrain Em F . We have concentrated here on the z = 4 UV LF as this is well defined over the widest range of 
luminosities. Improvements in the LF data, especially at the faint end, and at higher redshifts together with a better 
understanding of the astrophysics of galaxy formation, will allow us to place more stringent constraints on Em„ . 
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